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ABSTRACT 

Using N-body simulations we show that the process of formation of the brightest 
cluster galaxy through dissipationless galactic cannibalism can affect the inner cluster 
dark matter density profile. In particular, we use as realistic test case the dynamical 
evolution of the galaxy cluster C0337-2522 at redshift z = 0.59, hosting in its centre a 
group of five elliptical galaxies which are likely to be the progenitor of a central giant 
elliptical. After the formation of the brightest cluster galaxy, the inner cluster dark 
matter density profile is significantly flatter (logarithmic slope 0.49 </?< 0.90) than 
the original cusp (/? = I), as a consequence of dynamical friction heating of the massive 
galaxies against the diffuse cluster dark matter. In our simulations we have assumed 
that the cluster galaxies are made of stars only. We also show that the presence of 
galactic dark matter haloes can steepen the cluster central density profile. We conclude 
that galactic cannibalism could be a viable physical mechanism to reconcile - at least 
at the cluster scale - the flat dark matter haloes inferred observationally in some galaxy 
clusters with the steep haloes predicted by cosmological simulations. 
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1 INTRODUCTION 

Cold Dark Matter (CDM) cosmological simulations predict 
that the inner density profile of dark matter (DM) haloes is 
characterized by a cusp: Pdm(?~) oc r _/3 , with logarithmic 
slope P ~ 1 — 1.5 for r — * 0: the exact value of the slope and its 
universality are at the centre of a lively debate. For example, 
several authors (Navarro, Frenk & White 1996, NFW; Moore 
et al. 1998; Ghigna et al. 2000; Navarro et al. 2004) claim 
that the resulting profiles are universal, in the sense that /3 
is independent of the halo mass, while other authors suggest 
that the inner profile depends significantly on the slope of 
the power spectrum, so that the higher is the mass of the 
halo the steeper is the cusp (Subramanian, Cen & Ostriker 
2000; Ricotti 2003, and references therein). 

These predictions have been extensively tested by obser- 
vations both in galaxies (e.g., Salucci & Burkert 2000; van 
den Bosch et al. 2000; de Blok et al. 2001; Simon et al. 2003) 
and in clusters of galaxies (e.g., Smith et al. 2001; Sand, Treu 
& Ellis 2002; Kelson et al. 2002; Gavazzi et al. 2003; Lewis, 
Buote & Stocke 2003; Sand et al. 2004). The derived values 
of (3 range from ~ to ~ 1.3. The observational results - 
shallow slopes and intrinsic scatter - could represent a seri- 
ous challenge to the standard cosmological model, and exotic 
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scenarios have been also proposed to solve this problem (e.g. 
Spergel & Steinhardt 2000). 

In any case, the role of baryons must be better under- 
stood before we are forced to reject a successful paradigm. In 
fact, it is well known that on scales of the order of kiloparsecs 
baryons are important and their evolution may affect sub- 
stantially the DM distribution (see Binney 2004 for a discus- 
sion). For example, in the so-called "adiabatic contraction" 
approximation (Blumenthal et al. 1986; Gnedin et al. 2004), 
baryon dissipation effectively steepens the inner density pro- 
file of the host halo (Mo, Mao & White 1998; Kochanek & 
White 2001), thus exacerbating the contrast between theory 
and observations. 

Here, motivated by recent observational studies of 
galaxy clusters based on a joint lensing and dynamical anal- 
ysis - in which Sand et al. (2002, 2004) report inner loga- 
rithmic slopes in the range P = 0.38 — 0.99 - we focus on 
the cluster-sized DM haloes, which are considered to be less 



The uncertainties on f3 could be higher than those given by 
Sand et al. (2004) if the ellipticity of the cluster mass distribution 
is higher than estimated by the authors (cf. Dalai & Keeton 2004; 
Bartelmann & Meneghetti 2004). Additional observational work 
is currently under way to clarify those issues and measure the 
distribution of DM inner slopes for larger samples of clusters. 
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affected by baryon dissipation, with respect to galaxy-sized 
haloes. 

We note that the spiraling in of massive galaxies in clus- 
ters could flatten the cluster DM distribution, by heating it 
through dynamical friction (similar to the galaxy-globular 
cluster interaction at smaller scales; see Bertin, Liseikina & 
Pegoraro 2003). Recently, El-Zant et al. (2004), using N-body 
simulations, found that dynamical friction heating due to 
galaxies, modeled as rigid baryons clumps, is indeed effective 
in flattening the inner DM profile in clusters. Similar con- 
clusions, in a different context, were also reached by Ma & 
Boylan-Kolchin (2004), who explored with N-body simula- 
tions the effects of the dynamical evolution of (deformable) 
DM sub-haloes on the DM distribution of their host halo. On 
the other hand, the mass infall associated with the spiraling 
in of the galaxies towards the cluster centre deepens the cen- 
tral potential well. This dissipationless contraction, in con- 
trast with dynamical friction heating, can result in shrinking 
the cluster DM distribution and steepening its inner den- 
sity profile (see, e.g., Barnes & White 1984; Jesseit, Naab & 
Burkert 2002): the final density profile will be the result of 
the two competing effects (see also Section 4). 

Here we investigate the effects of dynamical friction on 
the cluster DM distribution in a realistic scenario, by sim- 
ulating the time evolution of the galaxy cluster C0337-2522 
at redshift z = 0.59 (ROSAT Deep Cluster Survey; Rosati 
et al. 1998). This cluster is characterized by the presence of 
five bright elliptical galaxies (Es) located within ~ 30 kpc 
from the centre, in projection. Nipoti et al. (2003b, here- 
after N03), using N-body simulations, showed that the five 
galaxies are very likely to merge and form a brightest cluster 
galaxy (BCG) by z — 0, in accordance with the predictions 
of the galactic cannibalism scenario (Ostriker & Tremaine 
1975; see also Merritt 1983; Dubinski 1998). In the present 
work we consider the evolution of the underlying cluster DM 
distribution, with the aid of additional higher resolution sim- 
ulations (see also Nipoti et al. 2003c). In our models the 
single galaxies, as well as the cluster DM, are represented 
by N-body distributions: this is particularly important, be- 
cause tidal stripping could remove mass from the individual 
galaxies, thus altering the infall process and its effects on the 
cluster DM halo. 



2 NUMERICAL SIMULATIONS 
2.1 Initial conditions 

Our simulations represent possible realizations of the dynam- 
ical evolution of C0337-2522 - a poor galaxy cluster, with 
five bright Es in its centre - from redshift z — 0.59 to the 
present. The observational data and the set up of the initial 
conditions are described in detail in N03. 

In the initial conditions of the simulations the cluster is 
represented as a live spherically symmetric DM distribution, 
hosting five identical massive galaxies modeled as spherical 
iV-body systems. For simplicity, we only consider simulations 
where each galaxy is represented by a single collisionless com- 
ponent (stars). This is justified under the assumption that 
the mass of the galaxies in the cluster centre is dominated 
by luminous matter (e.g. Treu & Koopmans 2004) , and that 
the extended DM halo is likely to be tidally stripped in the in- 
nermost region of the cluster (e.g. Natarajan, Kneib & Smail 
2002). 



Table 1. Parameters of the simulations. 
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First column: name of the simulation. r c y. initial cluster break 
radius (kpc). M c \: cluster mass (1O 13 M ). JVdm : total num- 
ber of DM particles. iVi um : total number of stellar particles. 
r DM f A*cl : best— fitting break radius of the final DM profile nor- 
malized to the initial break radius. /3: best— fitting logarithmic 
slope of the final DM profile, lcr uncertainties on ty>m,{ / r cl an d P 
are reported. Simulations 1.1, ls.l, 3.1, 17.1, and 20.1 are higher- 
resolution replicas of simulations 1, Is, 3, 17, and 20, respectively. 
Pairs of runs with the same cluster parameters (5a and 6, 18 and 
19, 20 and 21) differ in the initial centre-of-mass positions and ve- 
locities of the galaxies. In simulations Is and ls.l 50 small galaxies 
are added. 



2.2 Models 

For the cluster DM distribution we use a Hernquist (1990) 
model 



Pci(r) 



M cl r c 



27rr(r c i + r) 



(1) 



with Mci and r c i cluster total mass and break radius, re- 
spectively. Thus, p c i(r) ~ r' 1 for r — > 0, as in the NFW 
density profile, but, unlike the NFW profile, the total mass 
is finite, because p c i(r) ~ r" 4 for r — » oo. In the simula- 
tions the distribution function of the cluster has been cho- 
sen isotropic or radially anisotropic (see N03); r c i is either 
100 kpc or 300 kpc, and the total cluster mass is in the range 
4.8 x 10 13 - 3.06 x 10 14 M Q (see Table 1). 

At the beginning of each simulation the five galaxies are 
identical (isotropic or radially anisotropic) Hernquist mod- 
els: their stellar density distribution p» is described by equa- 
tion (1), with break radius r* ~ 2.2 kpc and total mass 
M* = 4 x 10 11 A#q. The initial centre-of-mass positions and 
velocities of the five galaxies are chosen by extraction from 
the cluster distribution function, constrained by the observa- 
tional data (see N03). Note that simulations with the same 
cluster parameters (namely, the pairs of simulations 5a and 
6, 18 and 19, 20 and 21; see Table 1) differ in the initial 
centre-of-mass positions and velocities of the galaxies. In sim- 
ulations Is and ls.l we consider an additional population 
of 50 smaller galaxies (each modeled as a Hernquist model 
with M, = 5x 1O 1O M and r* ~ 0.55 kpc) distributed in 
phase-space according to the cluster distribution function 
(see Fig. 1). 
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Figure 1. Simulation ls.l: snapshots of the initial (z = 0.59, 
left) and final (z = 0, right) projected distribution of the stellar 
particles in the inner cluster region. 



The total time of each simulation is ~ 6.1 Gyr, of the 
order of 100 Tdyn, where Tdyn i s the half-mass dynamical time 
of the initial galaxy models (we adopt fi m = 0.3, Qa = 0.7, 
and H = 65 km s" 1 Mpc" 1 ). 



2.3 Codes and numerical tests 

For most of the numerical N-body simulations (the first 12 
listed in Table 1) we used the Springel, Yoshida & White 
(2001) GADGET code, with cell-opening parameter a = 
0.02, minimum and the maximum time step A£ m i n = and 
At max = Tdyn/100, time step tolerance parameter a to i = 
0.05, and softening parameter e ~ 0.36r». Additional higher- 
resolution simulations (1.1, ls.l, 3.1, 17.1, 20.1) were ob- 
tained using the parallel code FVFPS (Fortran Version of 
a Fast Poisson Solver; Londrillo, Nipoti & Ciotti 2003). This 
code, based on algorithms introduced by Dehnen (2002), 
has been successfully tested and compared with GADGET 
(Nipoti, Londrillo, Ciotti 2003a; Nipoti 2003). In particular, 
we adopted the following values of the parameters: minimum 
value of the opening parameter # m i n = 0.5; softening pa- 
rameter e ~ O.lr,; initial time step At ~ Td yn /100. In all 
simulations stellar and DM particles have the same mass, 
and the total energy is conserved within 1%. 

Discreteness, two-body relaxation and softening are ex- 
pected to artificially alter the very inner DM density profile. 
To quantify numerical effects, we ran a few test simulations 
where the initial conditions are represented by the cluster 
DM halo only, for the adopted codes and values of the param- 
eters (softening, number of DM particles, time step, opening 
parameter, total elapsed time). The final profile of one of 
these test simulations (with the same parameters as simula- 
tion 1.1) is shown in Fig. 2 (empty symbols). In all the test 
simulations the Hernquist profile is preserved for r >0.05r c i. 
Thus, we only consider the properties of the cluster DM halo 
for r > 0.05r c i in the analysis of all the presented simulations. 
We note that in the initial conditions of the lowest-resolution 
of our simulations r ~ 0.05r c i contains ~ 270 particles. Thus, 
our choice of the minimum reliable radius is consistent with 
recent numerical convergence studies (Power et al. 2003, and 
references therein). 

Other possible numerical artifacts might be due to the 
small number of particles used to represent each galaxy (128 
to 2048, in the simulations run with GADGET). In partic- 
ular, it is likely that in our simulations we overestimate the 
galaxy mass that is tidally stripped (see Kazantzidis et al. 
2004) and this might affect the response of the cluster DM. 
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Figure 2. Final density distribution of the cluster dark, luminous 
and total matter in simulation 1.1. The dot-dashed curve is the 
initial (Hernquist) DM profile. The arrow indicates the minimum 
radius considered for the fit. Empty symbols refer to a test simu- 
lation run to exclude numerical artifacts (see Section 2.3). 



To explore this problem we ran simulations 1.1, ls.l, 3.1, 
17.1, and 20.1, which have the same initial conditions as sim- 
ulations 1, Is, 3, 17, and 20, respectively, but higher total 
number of particles (~ 10 6 DM particles, and 2048 to 8192 
stellar particles for each of the five galaxies) . We checked that 
all the relevant properties of the system evolve in the same 
way in the corresponding higher and lower resolution simula- 
tions. In particular the final cluster DM profile is practically 
indistinguishable in the two cases over the considered radial 
range r > 0.05r c i. 



3 RESULTS 

In all the simulations 3 to 5 out of the 5 massive galaxies 
merge before z = at the bottom of the cluster potential 
well, as a consequence of dynamical friction. Fig. 1 shows, 
as an example, the initial (z = 0.59, left) and final (z = 0, 
right) projected distributions of the stellar particles in the 
inner ~ 50 kpc for simulation ls.l. In any case, the merger 
remnant is similar in its main structural and dynamical prop- 
erties to a real BCG (see N03). Here we discuss the effects of 
this dissipationless multiple merging event on the underlying 
cluster DM distribution. 



3.1 Properties of the final dark matter density 
profile 

To quantify the combined effects of dynamical friction heat- 
ing and dissipationless contraction, we fit the final angle- 
averaged DM density profile using a generalization of the 
Hernquist profile, in the form (Dehnen 1993) 



Pdm(0 = 



^(r-DM+r) 4 -' 3 ' 



(2) 



where the inner slope f3 and the break radius rr>M are free 
parameters, and the reference density pdm,o is constrained 
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by the total DM mass. The initial cluster DM density profile 
(equation corresponds to the case /3 = 1, tdm = r c i, and 
Pdm.o = M c i/(27rr;?i). When computing the angle-averaged 
density profile of the final DM distribution, we determine the 
centre of the cluster using the iterative technique described 
by Power et al. (2003): the centre of mass of the system is 
computed recursively, considering particles within a sphere 
whose radius shrinks by 2.5 per cent at each step. We stop 
the iteration when the sphere contains ~ 1000 particles 2 . The 
best-fitting (3 and tdm for the final cluster DM distribution of 
each simulation are reported in Table 1; la uncertainties on 
the best-fitting parameters are calculated from A\ 2 = 2.30 
contours in the space /3 — roM. 

Fig. 2 shows the final dark, luminous and total matter 
distribution for simulation 1.1: the final cluster DM distribu- 
tion (dashed line) is centrally shallower than the initial Hern- 
quist profile (dot-dashed line). The modification is apparent 
at r < 20 kpc, where the contribution of the newly formed 
central galaxy (dotted line) to the total density (solid line) 
becomes relevant. Thus, the central cluster DM cusp after the 
formation of the BCG is flatter than the original pdm oc r^ 1 . 



1 




0.7 1 '''' 



0.01 0.1 

M |um /M lot 

Figure 3. Best-fitting logarithmic inner slope (top) and break 
radius r^M f (normalized to the initial break radius; bottom) of 
the final DM density distribution versus the ratio between the 
total mass in stars and total cluster mass (DM plus stars). Bars 
indicate lc uncertainties. 



We find the same qualitative behaviour of the final lu- 
minous and dark matter distributions in all the other sim- 
ulations, though simulation 1.1 represents one of the cases 
where flattening is more effective (the final best-fitting f3 
is ~ 0.5). Considering all of our simulations, we find best- 
fitting logarithmic slope f3 in the range 0.49 ^(3^ 0.90, with 
average = 0.71. Comparing the best-fitting final break 
radius rDM,f with the break radius r c \ of the corresponding 
initial Hernquist model, we find 0.74 < TDM.f /r c \ <0.95, with 
average (rDM,f/fci) = 0.85. We note that the best-fitting fi- 
nal break radius is not a direct measure of the concentration 
of the system, because it can be significantly affected by the 
change in the shape of the density distribution. In fact, the 
final radii containing 10, 50, 70 and 90 per cent of the total 
DM mass are always larger than 95 per cent of the corre- 
sponding initial values, even for quite small rDM,f/?"cl (e.g., 
simulations 1 and Is). Thus, in our models, deviations of 
roiii from r c i indicate a difference in the shape of the pro- 
file, as well as the value of (3 does. 

Simulation ls.l has the same initial conditions (cluster 
mass, radius and distribution function, and phase space co- 
ordinates of the centres of mass of the galaxies) as simula- 
tion 1.1, but in the first case 50 smaller galaxies are added to 
the five massive galaxies to represent the cluster population 
(see Fig. 1). We find that the main properties of the final 
matter distribution in simulation ls.l are not significantly 
different from those of simulation 1.1. We conclude that even 
a rather simple model, representing only the more massive 
galaxies in the initial conditions, is sufficient to capture the 
main physical processes. 

Thus, in the considered scenario we find that the most 
important effect of the formation of a BCG through galactic 
cannibalism on the underlying diffuse DM distribution is to 
produce a shallower final cluster DM distribution. In other 
words, dynamical friction heating is more effective than con- 
traction due to the mass infall. 



3.2 Physical interpretation 

In Fig. 3 we plot the final best-fitting (3 (top) and rDM,f/f c i 
(bottom) as a function of the ratio Mi um /M c i between the 
total mass in stars and the total cluster mass (DM plus 
stars). The general trend is that the final density profile is 
shallower when Mi um /M c i is higher. However, the distribu- 
tion of points in the plane Mi um /M c i-/3 (as well as in the 
planeMi um /M c i-rDM,f/?"ci) is characterized by a significant 
scatter: for example, simulations 1 and 3 have roughly the 
same M\ uul /M c \ ~ 0.04, but significantly different final inner 
slopes ((3 ~ 0.5 and f3 ~ 0.9, respectively). This suggests that 
the cluster mass is not a discriminant in determining the fi- 
nal density profile. We argue that the orbital properties of 
the five galaxies represent the main factor affecting the final 
slope of the DM cusp. This picture would be consistent with 
our findings, because simulations with the same or similar 
AZium/AZci (and cluster break radius) differ only in the ini- 
tial phase-space coordinates of the galaxies in the cluster. In 
general, galaxies on different orbits will heat the underlying 
DM at different radii and at different rates. Thus, it is not 
surprising to find a range of values of (3 in simulations with 
similar cluster parameters, but different initial positions and 
velocities of the galaxies. 

In order to quantify this effect, we present here the re- 
sults of a very simple exercise. We ran two additional test 
simulations, with the same number of particles, cluster pa- 
rameters and initial positions of the five galaxies as simula- 
tion 1, but different initial velocities. In one case each galaxy 
starts with half the initial kinetic energy it has in simula- 
tion 1 (and the same direction of the velocity vector), in the 
other the galaxies have null initial velocities. Interestingly, 
we find best-fitting final slopes f3 ~ 0.68 in the former case, 
and (3 ~ 1.29 in the latter, to be compared with (3 ~ 0.50 
for simulation 1. In other words, for fixed initial position of 

2 As a check, we also computed the centre of the system as the 
position of the lowest-potential energy particle, finding good agree- 
ment with the adopted method. 
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a galaxy, the smaller is its initial kinetic energy, the smaller 
the amount of energy it transfers to the cluster DM through 
dynamical friction heating. We also recall that the infall of 
the galaxies into the cluster centre deepens the cluster grav- 
itational potential well. If the amount of dynamical friction 
heating is quite small (as in the extreme case of the test simu- 
lation in which the galaxies start at rest), then dissipationless 
contraction dominates and the final cusp is steeper than the 
initial. 

We conclude that the orbital parameters of the galaxies 
are important in determining the final DM distribution, as 
clearly illustrated by the results of the two test simulations 
presented above. However, it must be stressed that these test 
simulations - intended only to isolate an important physical 
effect - are characterized by quite artificial initial conditions, 
which are not extracted from a distribution function and are 
not consistent with the observational phase-space constraints 
considered in our work (see N03). 

3.3 Comparison with observations 

How do the results of our simulations compare with ob- 
servationally determined cluster DM density profiles? The 
galaxy clusters observed by Sand et al. (2004) are suitable for 
this kind of exercise. These clusters (MS2137-23, Abell 383, 
RXJ1133, Abell 963, MACS1206, Abell 1201) are at redshifts 
lower than that of C0337-2522, they have total mass of the 
order of a few 10 14 M Q , and they host a dominant giant ellip- 
tical (with mass around 10 12 Mq) apparently relaxed at the 
bottom of the cluster potential well. Sand et al. (2004) found 
slopes 0.38 < ft < 0.99 with an average value of (ft) ~ 0.52 for 
the first three clusters, and an upper limit /3 < 0.57 for the 
others. It appears that our simulations are able to reproduce 
their observations. In particular, it is very suggestive to com- 
pare the final matter distribution of simulation 1.1 with the 
best-fitting model that Sand et al. found for MS2137-23 (cf. 
Fig. 2 with their Fig. 7) . The similarity of the DM and stel- 
lar distributions is striking: the best-fitting slope is /3 = 0.57 
for MS2137-23 and (3 = 0.50 for the end-product of simu- 
lation 1.1. Not only can the proposed mechanism reconcile 
the shallow observed slopes with the steep slopes obtained in 
DM only simulations, but also the various realizations of the 
cannibalism process could introduce additional scatter, thus 
helping to explain the observed scatter. 

Ultimately, when more observations are available and we 
can measure the distribution of the final slopes it might be 
possible to use the observed scatter to constrain the history 
of cannibalism in clusters. 



to the formation of a BCG. In our models, where the initial 
condition are realistic and the galaxies are deformable, we 
find final slopes in the range 0.49 < P < 0.90, with average 
(/3) = 0.71. The qualitative behaviour of our simulations is in 
agreement with what found by El-Zant et al. (2004) , though 
they obtained a quite small value of the final inner slope 
(/? ~ 0.35), which is likely to be due to the use of rigid 
clumps, combined with a high fraction of luminous (clumped) 
matter (Mi um /M c i ~ 0.06). Dynamical friction heating in a 
cluster with original inner density profile pdm k r" 1 can 
thus produce a flatter final inner slope for the cluster DM. 
Recent measures of the DM profiles of observed clusters can 
be interpreted taking into account this mechanism. 

As outlined in the Introduction, our results are based 
on the assumption that cluster galaxies are essentially made 
of stars. However, in a scenario in which the contribution of 
galactic DM haloes to galaxy masses is important, the in- 
spiraling galaxies might bring a significant amount of DM in 
the centre of the cluster. For instance, if we consider, as ex- 
treme case, that the massive galaxies in our simulations are 
made of DM only, the final DM cusp would be steeper than 
r , because the DM infall is dominant over dynamical fric- 
tion heating (cf. solid curve in Fig. 2). A similar scenario has 
been recently considered by Ma & Boylan-Kolchin (2004), 
who studied the interplay between a DM halo and its sub- 
haloes (see also Zhang et al. 2002). Ma & Boylan-Kolchin 
show with N-body simulations that the dynamical evolution 
of sub-haloes, depending on their mass and concentration, 
can either steepen or flatten the halo cusp. Thus, the deter- 
mination of the amount and distribution of DM in cluster 
galaxies is very relevant also to the problem of the cluster 
cusp. 

As a consequence, the discrepancies found between the 
observationally determined inner DM profiles and the predic- 
tions of CDM models must not necessarily represent a failure 
of CDM, at least at the cluster scale. Instead, they could re- 
flect our poor understanding of the formation and evolution 
of galaxies, where the physics of baryons is important. These 
processes have to be taken into account when predicting the 
properties of DM haloes of clusters, on scales of few tens of 
kiloparsecs or smaller. At the same time, our simulations - 
when interpreted in the context of structure formation (i.e., 
galaxies considered as sub-haloes of DM) - indicate that is 
quite remarkable to obtain a universal DM profile as a result 
of successive aggregations 3 , due to the requirement that dy- 
namical friction heating, potential well contraction, and mass 
increase all balance. This detailed balance should be further 
studied in cosmological simulations. 



4 SUMMARY AND CONCLUSIONS 

In this paper we explored, with the aid of N-body simu- 
lations, the effects on the inner cluster DM density profile 
of the formation of a BCG, through dissipationless multiple 
merging of pre-existing galaxies. The initial conditions are 
designed to reproduce the observed galaxy cluster C0337- 
2522 (z = 0.59). Although we focus on a single well observed 
system as a case study, we stress that the results of the sim- 
ulations could be generalized to poor galaxy clusters where 
the central giant elliptical was formed as a consequence of 
galactic cannibalism, with negligible dissipation. 

We found that the inner cluster DM profile is sensitive to 
the dynamical evolution of cluster galaxies, and in particular 
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